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RESEARCH MEMORAIilDUM 

SPATIAL BURNOUT IN WATER REACTORS WITH NONUNIFORM STARTUP 
DISTRIBUTIONS OF URANIUM AND BORON 
By Thomas A. Fox and Donald Bogart 


SUMMARY 

Spatial hiarnout calculations have been made of two types of water 
moderated cylindrical reactor using boron as a burnable poison to in- 
crease reactor life . Specific reactors studied were a version of the 
Submarine Advanced Reactor (SAE) and a supercritical water reactor 

(sew) , 

Burnout characteristics such as reactivity excursion, neutron- flux 
and heat-generation distributions, and uranium and boron distributions 
have been determined for core lives corresponding to a burnup of approx- 
imately 7 kilograms of fully enriched uranium. All reactivity calcula- 
tions have been based on the actual nonuniform distribution of absorbers 
existing during intervals of core life. Spatial burnout of uranium and 
boron and spatial build-up of fission products and equilibrium xenon 
have been- considered. Calculations were performed on the NACA nuclear 
reactor simulator using two-group diffusion theory. 

The following reactor burnout characteristics have been 
demonstrated : 


1 . A significantly lower excursion in reactivity during core life 
may be obtained by nonuniform rather than uniform startup distribution 
of uranium. Results for SOW with uranium distributed to provide con- 
stant radial heat generation and a core life corresponding to a uranium 
burnup of 7 kilograms indicated a maximum excursion in reactivity of 
2.5 percent. This compared to a maximum excursion of 4.2 percent ob- 
tained for the same core life when uranium was uniformly distributed 
at startup. Boron was incorporated uniformly in these cores at startup. 
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2. It is possible to approach constant radial heat generation dur- 
ing the life of a cylindrical core by means of startup nommiform radial 
and axial distributions of uraniim and boron. Results for SOW with non- 
uniform radial distribution of iiranium to provide constant radial heat 
generation at startup and with boron for longevity indicate relatively 
small departures from the initially constant radial heat generation dis- 
tribution during core life. Results for SAE with a sinusoidal distribu- 
tion rather than uniform axial distributions of boron indicate signifi- 
cant Improvements in axial heat generation distribution during the 
greater part of core life. 

3. Uranium investments for cylindrical reactors with nonuniform 
radial uranium distributions \diich provide constant radial heat gener- 
ation per imlt core volijme are somewhat higher than for reactors with 
imifonn uranium concentration at startup. On the other hand, uranium 
investments for reactors with axial boron distributions which approach 
constant axial heat generation are somewhat smaller than for reactors 
with uniform boron distributions at startup. 


UmODUCTION 

The use of distributed burnable poisons in power reactors to 
achieve increased core life and to eliminate a large part of shim con- 
trol has been imder consideration for some time. Appropriate distribu- 
tion of burnable poisons can permit fuel burnout to be balanced by poi- 
son burnout so as to maintain a reasonably small excursion in reactiv- 
ity. The advantages of using burnable poisons in reactors required to 
maintain specified total power outputs were pointed out in reference 1. 
The application of burnable poisons to naval reactors was illustrated 
in reference 2 and its possibilities for increasing reactor life and 
slu^ilifying control systems demonstrated. Reference 2 further suggested 
achievement of a more uniform heat flux and increased power densities by 
redistribution of the remaining core control rods present in a reactor 
such as the Submarine Thermal Reactor. 

The principle of distributed burnable poisons is being applied in 
the design of the Submarine Advanced Reactor (SAR) by Knolls Atomic 
Power Laboratory (KAPL) . In this connection MCA has cooperated with 
the Naval Reactors Branch of the Atomic Energy Commission and with KAPL 
in undertaking a study of achieving more uniform heat generation dis- 
tributions in high specific-power reactors incorporating burnable poi- 
sons. In addition, the variation of reactivity with operating time at 
a fixed reactor power as affected by local axial or radial spatial 
burnout has been determined. 
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Because the -use of distrihuted poison reduces the over-all amount 
of excess reactivity which is required to be controlled^ it is possible 
that core control rods may be entirely replaced by control rods at the 
reflector-core interface. Removal of control rods from the core makes 
possible the achievement of more uniform power distributions by nonuni- 
form startup distributions of uranium and poisons. 

The present report presents the resxilts of spatial burnout calcula- 
tions for two cylindrical water reactors incorporating natural boron for 
longevity purposes. Several variations in startup axial and radial dis- 
tributions of boron and uranium have been considered. The calculations 
have been performed on the NA.CA nuclear reactor simulator (ref. 3) using 
two-group diffusion theory. 


REACTOR LORGEVITI CRITERIA 

The use of btirnable poisons to increase reactor life affects reac- 
tivity excursion and heat generation distribution. Reactor design cri- 
teria regarding nonuniform startup distributions of uranium and boron 
are discussed. 

Reactivity exciarsion . - Because of the more rapid burnout of boron- 
10 relative to uranium- 235, a reactor goes supercritical after startup. 
However, the loss in reactivity due to burnout of urani\im-235 and accu- 
mulation of fission products resiilts in the reactivity passing through 
a maxim-um. This excursion in reactivity must be compensated for by in- 
sertion of control rods or by reduction in reflector effectiveness dur- 
ing core life. A reactor longevity criterion is, therefore, to keep 
this excursion in reactivity to a minimum for a given burnable poison 
and specified core life. 

The technique of "lumping" the burnable poison so that it shields 
itself at startup and makes available incremental quantities of poison 
as burnup proceeds has been considered in reference 4. It would seem 
possible by this means to flatten the reactivity excursion. In the 
present calculations, however, self- shielding effects have not been 
considered. What have been obtained are the burnout characteristics 
of a mmiber of initially critical reactors with different spatially 
distributed concentrations of urauium-235 and boron-10 at startup. 

The initial investments were chosen so as to achieve a uranium 
burnup of about 7 kilograms; this corresponds to a gross heat gener- 
ation of 133,000 megawatt hoirrs. 

Heat generation distribution . - Startup concentrations of uranium 
and boron are usioally designed to be uniformly distributed in the re- 
actor core. This would be conveniently obtained in practice by using 
identical irranium-boron elements equally spaced in the reactor core. 
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Heat generation distributions obtained for such reactors generally indi- 
cate peaks near the center of the core. These peaks ^ of course seri- 
ously reduce the total power available from a reactor operating at lim- 
iting fuel-element temperature. 

The objectionable peaks in heat generation can be removed by non- 
uniform distribuJ;ion of -uranium at startup. Techniques of initially 
achieving constant heat generation per xmit core volume by nonuniform 
distribution of fissionable material are discussed in references 5 and 
6. These techniques have been applied herein to indicate how a reactor 
designed for constant heat generation and incorporating a uniform con- 
centration of burnable poison at startup is affected by spatial burnout. 

Nonuniform startup distribution of burnable poison in the core can 
also improve the power distribution by altering relative spatial compe- 
tition with uranl-um for available neutron flux. However, the higher 
neutron fluxes near the center of the core cause more rapid b-urnout of 
boron here than in the remainder of the core-. The result is that during 
and at the end of core life, in reactors with initially uniform boron 
concentration, a significant quantity of boron remains in relative low 
fl-ux regions. This serves only to increase -uranium investment with no 
effect on reactor life. 

Therefore a reactor longevity criterion for nonuniform boron dis- 
tributions at startup is to produce (a) improvements in heat generation 
distribution and (b) savings in critical uranium investment as a result 
of smaller total quantities of boron necessary at startup to provide re- 
quired core life. 


SPECIFIC REACTORS 

Submarine Advanced Reactor (SAR) . - The first reactor considered 
corresponds to a version of SAR and consists of a cylindrical core 22 
inches in diameter and 44 Inches in height. It is cooled and moderated 
by -water and contains stainless steel - zirconium fuel elements . The 
axial and radial reflectors are mixtures of water and stainless steel. 
Reactor specifications and nuclear constants used in the SAR calcula- 
tions are given in table l(a) . Nominal reac-tor power is specified at 
65 mega-watts and average reactor temperature at 575° F| a 7 kilogram 
uranium bumup wo-uld therefore correspond to a core life of 2000 ho-urs. 

Supercritical Water Reactor (SCW) . - The second reactor considered 
is cooled and moderated by supercritical water. It consists of a cylin- 
drical core 30 inches in diameter and 30 inches in height and contains 
stainless steel fuel elements. The axial and radial reflectors are 
assumed to be effectively infinite thicknesses of water. Reactor speci- 
fications and nuclear constants used in the SCW calculations are given 
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In table l(b). Nominal reactor power is specified at 300 megawatts ; av- 
erage reactor core temperatin’e at 620° and average reflector tempera- 
ture at 480° F. A 7 kilogram uranium burnup for this reactor wo^ild 
therefore correspond to a core life of about 450 hours. 


REACTOR SPATIAL BURNOUT CALCULATIONS 

For cylindrical reactors constructed of fuel elements running the 
height of the core, improvements in axial heat generation distributions 
may be more practically achieved from a fabricatlonal viewpoint by a 
nonuniform axial distribution of b-urnable poison rather than uranium. 
However, in the radial direction, the possibility of using fuel elements 
with different uranium concentrations permits improvements in heat gen- 
eration distributions by nonuniform radial distribution of uranium. 

Cylindrical reactors inco3rporating the following distributions have 
been separately evaluated for the desired core life: 

1. Nonuniform concentrations of boron and uniform concentrations of 
uranium in the axial direction for SAR. 

2. Nonuniform concentrations of uranium and uniform concentrations 
of boron in the radial direction for SC¥. 

Spatial burnout characteristics such as reactivity excursion, core 
life, neutron flux, and heat generation distributions for the two water 
reactors have been calculated using two-group diffusion theory. The 
neutron diffusion equations, methods used for evaluating the two-group 
constants, and burnout relations are given in appendix B. (Symbols used 
are defined in appendix A.) 

In order to achieve the desired inranium burnup of about 7 kilo- 
grams, the burnout characteristics of a number of initially critical 
reactors with different spatially distributed concentrations of uraniiim- 
235 and boron-10 have been obtained. Approximate startup compositions 
for these reactors were determined by homogeneous hare-pile calciolations 
using assumed values of reflector savings (see appendix C). Reactivity 
exc\xrsion and core life for each reactor were determined by considering 
individual burnout periods which are successive intervals of the ex- 
pected life of the reactor core . The spatial flux distribution and 
burnout characteristics were assumed to remain constant during each in- 
terval of the reactor life. Relatively long burnup intervals may be 
considered without loss in accuracy because heat generation and fliix 
distributions change rather slowly with burnup as shown in appendix C. 
However, it is also shown in appendix C that longevity characteristics 
determined by homogeneous bare-pile calcilLatlons based on the variation 
with time of average rather than local core properties lead to serious 
inaccuracies . 



6 




MCA RM E55C16 


Calciilation Procedure 

The reactivity calculations have been performed on the MCA nuclear 
reactor sim-ulator. This electrical analog and the procedure for using 
it to solve reactor criticality problems have been described in refer- 
ences 3, 5, and 7. Although the reactor simulator employed is one- 
dimensional, solutions obtained include the effects of neutron leakage 
in the direction normal to the direction of solution. These effects 
are included by estimating reflector savings in this direction (see 
two-group equations in appendix B) . Solutions obtained on the simtila- 
tor automatically satisfy the required boundary conditions at the reac- 
tor center, core-reflector interface, and outer reflector boundary. 

Absolute flux levels, required only for equilibrium xenon determi- 
nation, were evaluated relative to the fluxes existing in the core at a 
position producing the nominal average power per unit core voliame. Av- 
erage equilibrium xenon concentration is approximately constant at these 
flux levels, so that the res-ults obtained are applicable to wide ranges 
of total core power and reactor lifetime. 

Both reactors considered are partly epithermal with 35 percent of 
the fissions for SAR and 20 percent of the fissions for SOW produced at 
energies above thermal at startup. Bpitheimal absorptions in -uranium 
and their effects on thermal fluxes and equilibrium xenon concentrations 
have been considered in the simulator spatial burnout calculations with- 
in the limits of the two-group approach as outlined in reference 8. 

Macroscopic cross sections and two-group parameters for a typical 
case of each reactor eire given in table II. The uranium-235 and boron- 
10 concentrations listed are averages over the spatial distribution ex- 
isting at the end of each individual burnup interval. Reactivity calcu- 
lations, of course, were based on the actual nonuniform distribution of 
absorbers existing during intervals of core life. Spatial burnout of 
uranium and boron and spatial buildup of fission products and equilib- 
rium xenon have been considered. Burnout relations used in computing 
spatial \xranl\mi and boron distributions are given in appendix B. 


REACTOR BURNOUT CHARACTERISTICS 

Spatial burnout characteristics of reactivity excursion, core life, 
and heat generation distribution determined for the axial direction for 
SAR and for the radial direction for SCW by the simulator spatial calcu- 
lations are presented. 
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Reactivity Variation 

The effective multiplication factor Kgpp as a function of uranium- 
235 hurnup is presented for both SAR and SOW in figure 1, (The circles 
indicate core b^umup intervals considered.) 

SAR . - Reactivity variations shovn in figure l(a) are for SAR with 
uniform axial distributions of boron and uranium at startup. Two cases 
with different initial investments of boron and uranium^ SAR I and 
bracket the prescribed core life of 7 kilograms uranium bumup. Larger 
reactivity excursions accompany increased core life. A reactivity ex- 
cursion of 2.5 percent was obtained for SAR which achieves closely 
the desired core life. This case required 30.0 kilograms of uranium-235 
and 202 grams of natural boron at startup. Starting compositions^ ex- 
cursion, and burnup for all cases are summarized in table III. Excess 
reactivity requirements have not been considered. 

Results in figure l(b) are for SAR cases III and IV with uniform 
distributions of \iranium and sinusoidal distributions of boron in the 
axial direction at startup. These startup boron distributions were 
chosen because they are proportional to the heat generation curve and 
so may be expected to approach xiniform minimal concentrations of boron- 
10 at the end of core life. It will be seen later that a startijp non- 
uniform boron distribution results in an in^jrovement in heat generation 
distribution by flattening the axial neutron fluxes. There is also a 
savings in critical -uranium investment because smaller quantities of 
boron are initially necessary to provide required core life. 

The 202 grams of natural boron of SAR II were distributed slniisol- 
dally in case SAR III. As shown in figure l(b), a core uranl-um burnup 
of 8.1 kilograms was achieved. To obtain the req^lired core life a 
fourth case, SAR IV, was computed in which the startup uranium invest- 
ment was reduced to 27.3 kilograms and Initial sinusoidally distributed 
boron req-ulrement was reduced to 148 grams. SAR IV resulted in a core 
burnup of 6.4 kilograms so that the required bumup of 7 kilograms is 
bracketed. The reactivity excin:*sion for the required burnup, interpo- 
lated from figtare l(b), is not affected by sinusoidal boron distribution 
in the axial direction. 

Results for SAR IV Indicate that required reactor life may be ob- 
tained by nommiform axial distribution of boron using less boron and 
less uranium and with about the same excursion in reactivity as in the 
uniformly distributed boron case. 

sew . - The reactivity variations sho-wn in figure l(c) are for SOW 
with uniform radial distributions of boron and uranium at startup. Two 
cases, sew I and II, were calculated. A reactivity excursion of 4.2 
percent was obtained for SGW II for a 7 kilogram uranium bumup. The 
sew core contains a lower concentration of uranium-235 than SAR so that 
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■burnup of a given mass of uranium represents a larger fraction of the 
xiranium concentration in SOW than in SAE. Inasmuch as the fractional 
horon-10 hurnup is proportional to the fractional uranlimi-235 hurnup 
raised to a power of 6.34 (see appendix B)^ horon-10 burns out more 
rapidly in SOW than in SAE and thus releases more excess reactivity 
after startup. SOW II required 25.1 kilograms of uranium- 235 and 251 
grams of natural boron at startup. 

Eestilts in figure l(d) are for SOW III in which boron is \miformly 
distributed and uranium-235 is distributed radially in nonuniform con- 
centrations so as to produce constant radial heat generation. SOW III 
required a critical mass of 29.5 kilograms of uranium and 217 grams of 
natiiral boron at startup. SOW III achieved a core burnup of 6-6 kilo- 
grams of \xranium-235 and a maximum reactivity excursion of 2.5 percent. 
The heat generation distribution for this case is practically uniform 
throughout core life as will be presently seen. 

Eesults for SC¥ III indicated that required reactor life may be ob- 
tained by nonuniform radial distribution of uranium using somewhat higher 
uranixim and lower boron investments and with a significantly lower ex- 
cursion in reactivity than in the uniformly distributed uranium case. 


Heat Generation Distribution 

Heat generation distributions for representative SAE and SOW cases 
at startup and at end of core life are presented in figure 2. The ra- 
tio of local heat release to average heat release per unit core volume 
h/h is shown as a function of reactor axial height for SAE and as a 
functioh of reactor radius for SOW. The distributions gradually assume 
final shape as burnout proceeds. 

SAE . - The distributions shown in figure 2(a) are for case SAE II 
and are representative of patterns obtained with unifoim concentrations 
of uranixmi and boron at startup. A maxlmum-to-minlmum heat release ra- 
tio of 2.8 at startup and 2.7 at end of core life are obtained. A sharp 
uptxirn in thermal flux at the core -reflector interface results in a 
large relative heat release ratio in the last few centimeters of core. 

In figure 2(b) are shown the distributions for case SAE IV with 
xmiform uranium and sinusoidal boron spatial concentrations at startup. 
These distributions are somewhat flatter than those obtained for the 
initially uniform boron cases. Maximum-to-minimum heat release ratios 
of 2.0 at startup and 2.5 at end of core life which are obtained in- 
dicate an improvement over SAE II. Most of the improvement occurs in 
the first half of core life. This Improvement gradually decreases as 
the net effect of burnout of boron and uranium no longer tends to flat- 
ten thermal neutron fluxes. 
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sew . - The distributions shown in figuire 2(c) are for case SOW II 
and are again representative of patterns obtained with uniform concen- 
trations of uranium and boron at startup. As in case SAR II, the dis- 
tribution improves slightly as uranium burnup proceeds preferentially 
in regions of power peaks. The reduced relative uranlm concentration 
causes reduced heat generation at these peak locations as burnout pro- 
ceeds. A maximum-to-minimum heat release ratio of 2.3 at startup and 
1.9 at end of core life are obtained. There is a sharp upturn in heat 
release ratio at the core -reflector interface region. 

The heat distributions shown in figure 2(d) are for case SOW III 
where uranium has been spatially distributed so as to produce constant 
radial heat generation at startup. The boron concentration was uniform 
and Identical to SOW I. The maximum-to-minimum heat release ratio is 
therefore 1.0 at startup. This ratio gradually changes to 1.2 by the 
end of core life. Thus, the radial heat generation distribution remains 
effectively constant throrighout core life and no sharp power gradients 
are present. 


Neutron Flux Distributions 

SAR . - Two-group axial neutron fl-ux distributions for case SAR IV 
operating at constant total power are shown in figure 3(a). The local 
fluxes relative to the average thermal flux at startup q 

presented as a function of axial height above the reactor midplane. Op- 
eration at a constant total power during core life causes the thermal 
fluxes to rise in the core as biornup proceeds. The fliixes for sinusoi- 
dal axial boron distributions aire somewhat flatter than the usml fluxes 
obtained with initial uniform concentrations of boron. The sharp upturn 
in thermal flux at the core-reflector interface results in a large power 
gradient here. 

sew . - Two-group radial neutron flux distributions for case SOW III 
operating at constant total power axe shown in figure 3(b). For this 
case a constant radial heat generation is obtained by nonuniform radial 
distribution of uranium. These flux distributions are representative. 


Uranium Distributions 

Spatial distributions of uranium-235 for individual bxirnup inter- 
vals considered are presented in figure 3. Local uranium concentration 

relative to average concentration at startup N^l^ are shown for the 

axial direction for SAR and for the radial direction for SCW. Birmout 
relations used are given in appendix B. 
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SAR . - Axial lorani-um distributions for cases SAR II and IV are 
shown in figures 4(a) and (b)^ respectively. The uniform distribution 
at startup gradually assumes the final nonuniform distribution. Aver- 
age values of uranium concentration over the core relative to startup 
concentration N^/nq are included for each burnout interval. 

sew . - Radial uraniimi distributions for cases SOW II and III are 
shown in figures 4(c) and (d)^ respectively. Distributions for SOW II 
in figure 4(c) are similar to those obtained for SAR. However^ uranium 
distribution necessary to produce uniform power concentration at startup 
for sew III;, indicated in figure 4(d) ^ retains its general shape during 
core life. 


Boron Distributions 

The corresponding spatial distributions of boron for the in- 
dividual bumup intervals for both SAR and SCW are shown in figure 
5. Local boron concentrations relative to average concentration at 
startup R®/Wq shown for the axial direction for SAR and for the 

radial direction for SCW. Burnout relations used are given in 
appendix B. 

SAR . - The axial boron distribution for case SAR II with uniform 
concentration of boron at startup is shown in figure 5(a). Konuniform 
spatial burnout of boron leads to the irregular distributions shown with 
significant amounts of boron remaining in the region of the core- 
reflector interface. 

The axial boron distribution for case SAR IV with sinusoidally lo- 
cated concentrations of boron at startup is shown in figure 5(b) . Boron 
burnout for this case proceeds at a much greater rate at the center of 
the reactor with the result that the initial sinusoidal distribution ap- 
proaches a relatively uniform distribution at end of core life. As in- 
dicated previously^ this results in more desirable power distributions 
and lower critical uranium requirements. 

SCW . - Radial boron distributions for cases SCW II and III are 
shown in figures 5(c) and (d) - Distributions obtained as burnout pro- 
ceeds are similar to those obtained with SAR with uniform boron distri- 
bution at startup. 

SAR V . - Results for the sinusoidal boron distributions for SAR 
suggested an inquiry into the extent to. which it is possible to approach 
constant axial heat generation per unit core volume at startup by non- 
uniform distribution of boron. It is recognized that a constant rate 
of heat release axially is not the most desirable from an engineering 
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viewpoint j since it voiild be possible to extract larger total amounts of 
heat from a reactor operating at constant fuel-element temperature in 
the axial direction. It is also recognized that although constant heat 
release at startup may be approached by nonuniform boron disposition^ 
the rapid btirnout of boron results in a gradual departure from constant 
heat release toward the end of core life. 

With these considerations in mind^ boron was distributed axially 
so as to obtain an initially flat heat generation curve. For this case 
SAR Vy the boron distribution necessary is shown in figure 6 with boron 
concentration effectively zero at the core -reflector interface. Qlie re- 
sulting heat generation curve, also shown in figure 6, is flat except at 
the core -reflector interface. The critical uranium requirement for case 
^ SAR V is 55 kilograms and natural boron requirement is 620 grams, both 

3 considerably greater than req'uirements for cases SAR I to IV. The life 

2 of this core, which would be considerably greater than a burnup of 7 

J kilograms of uraniimi, has not been evalimted. The heat generation at 

i the interface rises significantly because of the feed-in of thermal 

neutrons from the reflector. An effective boron curtain at the inter- 
face woiild be required to suppress thermal fluxes here which in turn 
would require greater, amounts of uranium to maintain criticality. How- 
ever, improvements in heat generation are possible with boron distribu- 
tions intermediate between that of cases SAR IV and V. 


CONCLUSIONS 

From two-group spatial burnout calculations of two types of water 
moderated cylindrical reactor incorporating boron and uranium concen- 
trations in uniform or nonuniform spatial distribution at startup, the 
following reactor burnout characteristics have been demonstrated for a 
given core life: 

1. A significantly lower excursion in reactivity during core life 
may be obtained by nonuniform rather than uniform startup distribution 
of uranium. Results for a supercritical water reactor with uranium 
distributed to provide constant radial heat generation and a core life 
corresponding to a uranlimi b'urnup of 7 kilograms indicated a maximum 
excursion in reactivity of 2.5 percent. This compared to a maxlmm 
excursion of 4.2 percent obtained for the same core life Tdien uraniima 
was uniformly distributed at startup. 

2. It is possible to approach constant radial heat generation dur- 
ing the life of a cylindrical core by means of startup nonuniform radial 
and axial distributions of uranium and boron. Results for a supercriti- 
cal water reactor with nonuniform radial distribution of ^iranium to pro- 
vide constant radial heat generation at startup and with boron for lon- 
gevity indicate relatively small departures from the initially constant 
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radial heat generation distribution during core life. Results for a 
version of the Submarine Advanced Reactor with sinusoidal rather than 
uniform axial distributions of boron indicate significant improvements 
in axial heat generation distribution during the greater part of core 
life. 


3. Uranium investments for cylindrical reactors with nonuniform 
radial uranitim distributions which provide constant radial heat gener- 
ation are somewhat higher than for reactors with uniform uranium con- 
centration at startup. On the other handj uranium investments for re- 
actors with axial boron distributions which approach constant axial heat 
generation are somewhat smaller than for reactors with imiform boron 
concentration at startup. Considerable design latitude would seem to 
be available for meeting reactor heat generation and core life require- 
ments through the lose of nonuniform startup distributions of uranium 
and boron with approximately the same uranixmi Investments as for reac- 
tors in \dilch uranium and boron are uniformly distributed at startup. 


Lewis Plight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, March 23, 1955 
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APPEOroiX A 
SYMBOLS 

The following symbols Are used in this report: 
buckling constant 

K neutron multiplication constant 

effective mioltiplication factor 

p 

Lf mean-squEire slowing-down length for fast neutrons 

2 

mean- square thermal diffusion length 
N atom concentration per unit volume 

P total fissions per unit volume of reactor core 

p^j^ fraction of neutrons absorbed thermally 
t time 

V neutron velocity 

a ratio of radiative capture cross section to fission cross section 

P fraction of total fissions/(cc) (sec) occurring thermally 

macroscopic transport mean free path 

V neutrons produced per fission 

2 macroscopic neutron cross section 

cj microscopic neutron cross section 

(p neutron flux 

Subscripts : 

A absorption process 

F fission process 

f fast neutron group 

th thermal neutron group 
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APPENDIX B 


DIPFUSION EQUATIONS AND BURNOUT RELATIONS 

Two-group diffusion egmtlons . - The two-group neutron-diffusion 
eqiiations for the cylindrical reactor configiorations considered are: 


Core- 




V^'Pth 


K. 


1 'P-P 

■(1 - Pth) - pj ^ 


’’f . *-IR,th '^h 




-“th 


'Pth + 


X 


TR,f ^th 


"th 


^TR;,th l| 


«Pf - B^cpth = 0 


V^cPf - i spj. - B^qif = 0 

4 


Reflector < 


V^tp- 


th - ‘Pth + ^f - B^^th = 0 

^TR,th l| 


Axial solutions to these equations have been obtained for SAR and 
radial solutions have been obtained for SC¥. Epithermal fissions are 
included and are considered to occur in the fast group. Leakage of 
neutrons in the direction normal to the direction of solution is also 
included by consideration of the buckling constant for the equivalent 
unreflected dimension in the normal direction. Solutions obtained on 
the electrical-analog nuclear reactor simulator automatically satisfy 
the usual boundary conditions at the core center, core-reflector in- 
terface, and outer reflector boundary. 

The procedure for evaluating nuclear constants for use in the two- 
group equations is described in references 9 and 10. The procedure is 
briefly outlined as follows : 

The fast group constants K^, ^ for core and reflector 

have been averaged from fission to thermal energies by the flux dis- 
tribution in an Infinite medium of the same composition according to 
the Femnl age slowing-down model. Dependence of this flux distribu- 
tion on the fission spectrijm is included. 

p 

For water, L^ is based on the experimental value of 33 square 

centimeters at room temperature (p = 1 gm/cc) and is taken as in- 
versely proportional to the square of the water density at higher 
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p 

water temperatiires . Values of for the core and reflector were com- 

puted by the following expression given by KAPL for hydrogenous media; 



(l + V)^ 

(p + xV)2 


where x is an equivalence factor which depends upon the particular 
metal used in conjunction with water, V is the metal-water volume ra- 
tio, and p is the density of water. The equivalence factors used 

were x^^ = 0.40, and = 0.38. In cores where two metals were 

present, a weighted average was obtained for x based on the volumetric 
concentration of the metals. 

2 

The thermal group constants have been evaluated 

by weighting local values according to the Maxwellian distribution of 
neutron flux for the temperature. The local absorption cross sections 
of all reactor materials with the exception of xenon- 135 are assumed to 
be inversely proportional to neutron velocity v. 

By means of the method of reference 11, the effects of chemical 
binding of hydrogen in water have been considered in evaluating 

Experimental values of ffg as measured for water are used to calciilate 
local values of cj^. The quantity for the mediimi is then 

evaluated by weighting local values by the Maxwellian distribution of 
neutron flux. 


Uranium- 235 burnout . - Uranium bijrnup rate is given by 

.U 


dt ^A,th‘Pth + 




T.th 


°"p,th 


= pP(l + o^h) + (1 - + °^f) 


or burnup is 




p(l + O^th) + (1 - P)(l + c^f)]pt 


In the present calculations, the core life was specified as that 
corresponding to a uranium burnup of 7 kilograms of uranlum-235. As 
such, the final value of Pt corresponds to 133,000 megawatt hoirrs. 

The assumption was made that the ratios of radiative capture to fission 
in uranium-235 for thermal and fast groups are equal, that is, 
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Boron-10 burnout . - Boron-10 burnup rate is given by 


dN® yB 
“ % 


th'Pth + 



with 


'Pth - 


PP 


"•F^th 

^ pp „B Cl - p)p 
" dt %;,th _u „u 

-^Fjth ^F,f 


or 




rB 




B 


B 


p + (1 - p) 


U 


L'^F.th '^F,f 


nU 


dt 


But 


Therefore 




0? 


p(l + + (1 - P)(l + a^) 


0? 


Ft 


^ P(1 + ^ (1 - p)(l + a^) 


N; 


B 


B 


B 


0 


K 


IjU 

0 


^A,th 


P(1 + *Hh) + (P " P)(P + “-f) 


In the present calculations, the assumption was made that the ra- 
tios of microscopic absorption cross sections for boron- 10 and uranium- 
235 for thermal and fast groups are equal. With the further assumption 
that the exponent above reduces to th/'^A th “ 6 * 34:. 

Fission products other than xenon . - The fission products other 
than xenon have been combined. An average fission product with a 
thermal absorption cross section 100 barns at 0.025 ev is 

assumed to be generated per fuel atom fissioned. Its cross section 


BH i anar 
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is assumed to follov the l/v law. Burnup of these fission products with 
reactor operating time is neglected. Fission product generation rate is 
given hy 


dM^ 

dt 


~ ^F,th‘*’th 


+ = PP + (1 - 3)P 



= Pt 


Equilibrium xenon. - Equilibrium concentration of xenon~135 is 


given by 




/ep 




^Xe 




°A^th^th ^Xe ^?thP^ 




F,th 


where 

yXe yield per fission = 0.067 

average over Maxwellian for thermal group 
decay constant = 2.1«10“^ per second 


Epithermal burnup . - In order to evaluate two-group epithermal 
burnup of uranium-235 and boron- 10 Independently, the average values 

of and cr^ f/°A f burnup expressions must be 

known for the reactors considered. The assumptions that a^. = 

and °A f/°A f = th/°A th present calculations, were 

qualitatively checked using the Goertzel-Selengut model. Inasmuch as 
the Goertzel-Selengut model is a more appropriate slowing-down model 
for hydrogenous media, values of these fast parameters were obtained 
by this model for an infinite mediijm closely corresponding to case 
SAR I startup composition; average values over the epithermal region 

of df = 1.85 and' = 0.80 °A^th/'^A,th 'were calculated. 

Use of these values in the burnup expressions indicates that 
xiranium-235 burns out 4.5 percent faster than for the entirely thermal 
assumption and therefore boron-10 also burns out faster (for case SAR I 
at 2000 hours, boron-10 burnup is increased by 5.0 percent). The net 
effect on reactivity for case SAR I is estimated to make the reactor 
considering epithermal effects only 0.1 percent more reactive than the 
reactor calculated under the foregoing thermal assumptions . 
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APPENDIX C 


HOMOGENEOUS BAEE-PILE CALCULATIONS 

The two-group criticality equation for a homogeneous bare reactor 
is given by 


^h^th ^ - Pth) 

(l + l|b2)(i + l|j^b 2) (l + l|b2) 


where for cylindrical reactors 


B 


2 




and 


Eg = Ec + reflector radial savings 
Hg = Hq + 2 (reflector axial savings) 

Criticality calculations have been made using this equation and 
the radial and axial reflector savings specified in table I. These 
calculations served as a guide for the spatial burnout solutions in 
all cases except SCW III with uniform radial heat generation. For 
baxe-pile calculations^ burnouts were based on the fluxes in the av- 
erage -unit volume of the reactor core. Coir^parison of the variation 
of reactivity with iiranium burnup for SAE I and II as calculated by 
bare pile and spatial burnup techniques is shown in fig\xre 7. Hiere 
is a significant difference between the values of K^ff for the two 
calculations . 

A large part of the difference between the values of for 

the spatial burnout cases (calculated on the simulator) and the homo- 
geneous bare -pile calculations (based on the fluxes in the average unit 
voliime of the reactor core) shown in figure 7 is due to the procedure 
used in calculating the boron- 10 burnup. Because boron-10 burnup is 
related to uraniutm-235 burnup raised to a power of 6.34^ the actual 
boron -10 burnup is not the average uranium burnup raised to this power ^ 
but the sum of the local \iranium burnups, each raised to this power. 

The last procedure always results in a smaller boron burnup or more 
boron left. Therefore, reactors calculated on the simulator are less 
reactive than reactors calculated by homogeneous bare-pile procedure. 
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To illustrate these effects, the average boron-10 concentration at 

"D 

a core burnup of 6.64 kilograms, designated has been evaluated by 

the two methods. A comparison of as calculated on the simulator 

and by homogeneous bare -pile considerations is given in the following 
table : 


Method of calculation 

t 

1 atoms \ 
\cc core/ 

^eff 

(at time t ) 


SAR I 

SAR II 

SAR I 

SAR II 

Simulator 

1.50X1018 

1.93x1018 

0.988 

1.007 

Homogeneous bare pile 
Homogeneous bare pile with 

1.29 

1.70 

1.000 

1.015 

simulator ^ 

1.50 

1.93 

.995 

1.011 


The results in the table indicate that the difference in method of 
calculation of boron-10 accounts for about half of the difference in re- 
sultant values of The remaining difference between the simulator 

and bare-pile values of Keff probably due to the actual spatial 
variation of the thermal absorbers considered in the simulator solution. 
It should be noted that the homogeneous bare-pile calculations give good 
values of for reactor cores having uniform concentrations of 

strong thermal absorbers. However, values of as determined by 

the two types of calculations differ considerably with more extreme 
variations in spatial concentrations of thermal absorbers. For example, 
in sew III with large radial variation in uranium concentration, a homo- 
geneous bare-pile calculation has little meaning. 

As previously mentioned, the fltix distribution at the start of each 
burnout interval was assumed to hold for the entire interval. To check 
the sensitivity of this assumption, the reactivity after a uranium 
burnup of 6.64 kilograms was evaluated on the simulator for SAR I and 
II applying the flux distribution at reactor startup for the entire 
burnup period. These values are separately indicated in figure 7 at 
a core burnup of 6.64 kilograms and agree very closely with the smaller 
burnup intervals. These results indicate that the starting nontiniform 
spatial flux distributions are little affected by burnout and that 
rather large burnout Intervals may be considered with small loss of 
precision. 
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TABLE I. - EEACTQR SPECIFICATIONS AND CONSTAIJTS 


(a) The Submarine Advanced Reactor (SAE) 


T = 5750 Pj kT = 0.0495 ev: u,. = 19.12 
av ’ th ■ 

p^2° = 0.731; = 6.40; pSS-304 ^ ^ 


Con 5 )Osition by volume 

Core 

Axial reflector 

Radial reflector 


SS-304 

( 22 - in. dlam.; 

(lO-in. thickness; 

(l4-in 

thickness ; 



44-in. height) 

J 

)-cm reflector saving) 

12-cm reflector saving) 



HgO 0.600 


0.900 

^2° 

0.350 

Cr 0.202 

SS-304 .148 

SS-304 

.100 

SS-304 .650 

Ni .085 

Zr .244 






Fe . 713 

U .008 










1.000 


1.000 


1.000 


Atomic concentration 

. N, atoms/cc 




Atom 

Core 

Axial reflector 

H 


0.02936x10 

24 


0.04403X102^ 

0 


.01468 



.02201 

Zr 


.01034 






Cr 


.00258 



.00174 

Ni 


.00109 



.000734 

Fe 


.00911 



.00616 


Thermal cross sections (u,, = 

' th 

19.12) 




Atom 

Core 

Axial reflector 



°S' 


^TR’ 

2a ^ 


Ztr^ 


barns 

barns 

l/cm 

l/cm 

l/cm 


l/cm 

H 

.2246 

^23.3 ( 0 ^) 

0.00584 

0.6841 

0.00877 

1.0259 

0 

.0007 

4.0 

0 

.0563 

0 


.0844 

Zr 

.2457 

9.0 

.00225 

.0920 




Cr 

2.037 

4.2 

.00466 

.0107 

.00314 

.0072 

N1 

3.151 

17.3 

.00303 

.0186 

.00205 

.0125 

Fe 

1.678 

10.9 

.01354 

.0978 

.00915 

.0661 

^235 

446.8 

8.2 

(t) 

Negli- 





380 



gible 




rIO 

F 

2830 

4.0 

(b) 

Negli- 








gible 




Fission product 


0 

(b) 

0 





^Averaged over Jfexwelllan including chemical binding effects . 
^Varies with burnout . 




22 


mCA BM E55C16 


TABLE I. - Concluded. REACTOR SPECIFICATIONS AND CONSTANTS 


Core: 

Reflector : 


(B) The Supercritical Water Reactor (SCW) 

T^^ = 620° F; kT = 0.052 ev; u^j^ = 19.07; P®2° = 0.71; 
T^^ = 480° F; kT = 0.045 ev; u^^^ = 19.22; P^2° = o.83. 


8.03 


Composition by volume 

Core 


Axial reflector 

Radial reflector 

SS-347 

(30- in 

. dlam; 

(Effectively infinite 

(Effectively infinite 


30-in. 

height) 

thickness j 7 

.5 -cm 

thickness; 7. 

5 -cm 




reflector saving) 

reflector saving) 


HgO 

0.884 


HgO 1.000 

HgO 1.000 

Fe 0.690 

SS-347 .116 







Ni .085 

U 

0 







Cr .195 


1.000 







Mn .021 









Nb .009 









1.000 




Atcmic concentration 

, N, atoms/cc 



Atom 

Core 

Reflectors 


H 


9A 

0 ,04220 kio 


0 

05552 >d0^^ 


0 


.02110 



.02776 



Fe 


.006854 






Ni 


.0009057 






Cr 


« .001840 






Mn 


.0001933 






Nb 


.00005716 





Core thermal cross sections = 19.07) 


Atom 


^k> 

®s^ 


Ik, 

Ztr^ 




barns 

barns 


l/cm 

1/cm 

H 



0.2195 

®'22..1 

0.00926 

0.935 

0 



.0007 

4.0 

Negligible 

.0808 

Fe 



1.684 

10.9 

.01154 

.0740 

Ni 



3.119 

17.3 

.00283 

.0154 

Cr 



1.985 

4.2 

.00370 

.0076 

Mn 



8.73 

2.3 

.00169 

.0004 

Nb 



.762 

5.0 

.00004 

.0003 

u235 



440.0 

8.2 

(b) 

Negligible 

bIO 



2760.0 

4.0 

(b) 

Negligible 

Fission product 

69.3 

0 

(^>) 

0 

1 


^Averaged over Maxwellian including chemical binding effects. 
Waries with burnup. 




TABLE II. - MACROSCOPIC CROSS SECTIONS AND PARAMETERS' 


Parameter 

SAR II 

SOW II 


1 Burnupj kg U235 

Axial 



Burnup , 

kg u235 



Radial 


0 

1.66 

3.32 

4.98 

6.64 

reflector 

0 

1.5 

3.0 

4.5 

6.0 

7.0 

reflector 


3.504 

3.504 

3.504 

3.504 

3.504 

3.663 

3.899 

3.899 

3.899 

3.899 

3.899 

3.899 

4.13 

t2 2 

Lf, cm 

92.4 

92.4 

92.4 

92.4 

92.4 

69.5 

67.5 

67.5 

67.5 

67.5 

67.5 

67.5 

47.9 

^TR,th’ 

1.042 

1.042 

1.042 

1.042 

1.042 

.836 

.9180 

.9180 

.9180 

.9180 

.9180 

.9180 

.4747 

^th 

.676 

1.394 

.694 

1.416 

.712 

1.430 

.727 

1.428 

.740 

1.415 

.950 

.792 

1.300 

.808 

1.336 

.821 

1.354 

.832 

1.355 

.842 

1.343 

.848 

1.324 

.950 



Kf 

HrtO+SS , 

1.403 

1.422 

1.432 

1.430 

1.418 

1.267 

1.297 

1.312 

1.314 

1.303 

1.287 



Z ^ , l/cm 

A,th ' ^ 

.02930 

.02930 

.02930 

.02930 

.02930 

.02310 

.02609 

.02609 

.02609 

.02609 

.02609 

.02609 

.01208 

nlO 

%tb’ 

..235 

, l/cm 

.02088 

.01495 

.01016 

.00702 

.00484 


.02000 

.01371 

.00927 

.00618 

.00406 

.00303 


.11083 

.10469 

.09855 

.09242 

.08628 


.07382 

.06943 

.06504 

.06065 

.05626 

.05307 


A,th 















0 

.000831 

.001662 

.002492 

.003323 


0 

.00051 

.00102 

.00153 

.00204 

.00238 


^Xe ^2^ l/cm 

A,th' A,th ' 

.0727 

.0713 

.0703 

.0696 

.0689 


.0650 

.06524 

.05966 

.05870 

.05866 

.05879 


2^® , l/cm 

.00806 

.00746 

.00693 

.00643 

.00595 


.00480 

.00453 

.00388 

.00356 

.00330 

.00312 


A,th 














-jr235 p 

N *10^ ^ atcaas/cc 

2.800 

2.645 

2 .490 

2.335 

2.180' 


1.890 

1.780 

1.665 

1.555 

1.440 

1.360 


■olO 














B 20 , 

N 'lO , atoms/co 

.0833 

.0597 

.0406 

.0280 

.0193 


.0818 

.0561 

.0379 

.0253 

.0166 

.0124 


*^eff 

1.000 

1.019 

1.025 

1.020 

1.007 


1.000 

1.026 

1.040 

1.040 

1.022 

1.001 



Averaged over the nonuniform spatial distribution of absorbers existing at period of core life. 


CO 

OQ 
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TABLE III. - INITIAL INVESTMENTS, MAXIMUM REACTIVITY EXCURSION AND 
CORE LIFE FOR REACTOR CORES CONSIDERED 


Case 

Fully enriched 

y 235 startup 

requirement , 
hg 

(I) 

Natural horon 
startup requirement, 
g 

Maximum excursion 
reactivity, 
percent 

Reactor core 
life hurnup 
of u235^ 
kg 

SAR I 

28.2 

173 

1.65 

5.8 

SAR II 

30.0 

202 

2.50 

7,2 

SAR III 

29.5 

202 

2.70 

8.1 

SAR TV 

27.3 

148 

1.85 

6.3 

SAR V 

55.0 

620 

— 

— 

SC¥ I 

24.0 

217 

3,10 

5.6 

sew II 

25.1 

251 

4.18 

7.0 

sew III 

29.5 

217 

2,55 

6.6 


^'Exclusive of excess reactivity requirements. 
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0 2 4 6 8 10 

Core burnup, kg 

(b) Cases SAE III and IV. 

- Reactivity variation with core uranium burnup. 




Figure 1. 




Effective multiplication factor 
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Core turnup, kg 
(d) Case SCW III. 


Figure 1. - Concluded. Reactivity variation vith 
core uranium turnup. 
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0 10 20 30 40 50 60 

Reactor axial height, z, cm 

(h) Case SAR IV. 

Figure 2. - Heat generation distributions. 




Local to average heat releas 
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(c) Case sew II. 
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0 10 .30 40 

Reactor radius, r, cm 

(d) Case SCW III. 

Figure 2. - Concluded. Heat generation distributions. 
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0 10 20 30 40 50 60 70 80 

Reactor axial height, z, cm 

(a) Case SAR IV. 

Figure 3. - Relative neutron flux distributions. 





Relative neutron flux, <p/p 
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.1 I I I I L-J I 1 I I 

‘ 0 10 20 30 40 50 60 70 80 

Reactor radius, r, cm 

(■b) Case SCW III. 


Figure 3. - Concluded. Relative neutron flux distributions. 
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0 10 20 30 40 50 60 

Reactor axial height, z, cm 

(h) Case SAR IV. 

Figure 4. - Urani\am-235 distributions. 






Local to initial average uranium concentration 
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Reactor radius, r, cm 


(d) Case SCW III. 


Figure 4. - Concluded. Uranium- 235 distributions. 
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Figure 5. - Boron- 10 distributions. 
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0 10 20 30 40 

Reactor radius, r, cm 


(d) Case SCW III. 


Figure 5. - Concluded. Boron-10 distributions. 





Figure 6. - Startup boron distribution for initially uniform heat genera- 
tion. SAR V- 
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Effective multiplication factor, 





Figure 7. - Criticality variation with 


core uranium burnup. 
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